Electric-fleld-driven polymer entry into asymmetric nanoscale channels 
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The electric-field-driven entry process of flexible charged polymers such as single stranded DNA 
(ssDNA) into asymmetric nanoscale channels such as a-hemolysin protein channel is studied the- 
oretically and using molecular dynamics simulations. Dependence of the height of the free-energy 
barrier on the polymer length, the strength of the applied electric field and the channel entrance 
geometry is investigated. It is shown that the squeezing effect of the driving field on the polymer 
and the lateral confinement of the polymer before its entry to the channel crucially affect the barrier 
height and its dependence on the system parameters. The attempt frequency of the polymer for 
passing the channel is also discussed. Our theoretical and simulation results support each other and 
describe related data sets of polymer translocation experiments through the a-hemolysin protein 
channel reasonably well. 

PACS numbers: 87.15.A-, 36.20.Ey, 87.15.H- 



I. INTRODUCTION 

Polymer translocation - passage of a polymer through 
a pore in a membrane - is a very ubiquitous and vital 
process in biological environments. In the past decade 
it has also been used as a tool for single-molecule stud- 
ies [l|. Better understanding of related biological pro- 
cesses, DNA sequencing [2'], examining the existing the- 
ories on the statics and dynamics of confined polymers 
[1, 3 , direct study of nucleic acid-protein interactions Q 
and nucleic acids secondary structure f6s] are examples of 
present and potential applications. As a polymer passes 
through a pore, the number of its possible conformations 
and hence its entropy decreases. To overcome this en- 
tropic barrier, a voltage difference is often applied in the 
case of charged polymers. During the polymer passage, 
the ion current through the channel is blocked showing 
that the polymer is inside the channel 0|. Time length 
of polymer passage through the pore [8|, |9| , ion current 
during the polymer passage [l^, [HI and the time inter- 
val between consecutive passages (entry time) [Tll - [29| are 
three measurable quantities in polymer translocation ex- 
periments. These quantities contain information about 
the polymer and its interaction with the channel (30| . 
Despite the two first quantities, the polymer entry time 
into the channel is less investigated. 

The entry process of the polymer into the channel may 
be diffusion- or barrier-limited. It has been observed 
in the experiments that the entry time as a function of 
the inverse of the applied voltage is linear (exponential) 
in diffusion-limited (barrier-limited) regime hj. Fluo- 
rescence microscopy has revealed that in the diffusion- 
limited regime, there is a region of very low polymer 
concentration around the channel entrance. In this re- 
gion, the electric field drift dominates over the diffusion 



and the polymer is carried to the pore in a very short 
time. This is continued with immediate passage of the 
polymer through the pore p^ . In this regime, diffusion- 
convection (Smoluchowski) equation can explain the be- 
havior of the polymer outside the channel and its flux 
through the pore very well 14 17J. 

In the barrier-limited regime, the entry process of the 
polymer into the channel is very time consuming. The 
free-energy barrier against the polymer entry into the 
channel [l^, E^, [13 and the polymer attempt frequency 
for passing this barrier (l6l . [2l| are believed as two main 
quantities determining the entry time [31 . It has been 
shown that in this regime the entry time depends on the 
polymer length in contrast to the diffusion-limited regime 

Positional distribution of charges on the channel sur- 
face, electro-osmotic flow of counterions, interaction of 
the charges on the polymer with the channel medium 
which is of low dielectric constant and compression of the 
cloud of the polymer counterions in the course of poly- 
mer passage through the channel have been mentioned 
as other possible relevant factors [ll|, HI, . 

Enhancement of the frequency of polymer entry is of 
great interest in polymer translocation experiments. It 
provides the possibility of polyrner translocation studies 



in lower salt concentrations [2J|, lower applied voltages 
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[l^ and smaller sample amounts [T^. For this purpose, 
numerous procedures such as applying a salt gradient 
across the channel [l^, HI], charge manipulation of the 
protein channels by site-directed mutagenesis [l^ and 
layer deposition on solid state channels [11] have been 
used. 

In this paper, we investigate the polymer entry time 
in the barrier-limited regime for two different geometries 
of the channel entrance. A simple scaling theory as well 
as the results of coarse-grained molecular dynamics (MD) 
simulations are presented. We show that the electric field 
compresses the polymer on the wall before its entry to the 
channel. We calculate the free-energy of the polymer in 
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FIG. 1: Schematic of an asymmetric channel consisting of two 
coaxial cylinders of different diameters and a polymer which 
is released from equilibrium in front of the channel (upper 
panel). Schematic of the free-energy landscape of the polymer 
when it enters from the cis side of the channel without and 
with applied driving electric field (lower panel). 



the compressed state and then find the dependence of the 
height of the free-energy barrier on the polymer length 
and strength of the electric field for the two geometries of 
the channel entrance. Results of the theory and the sim- 
ulation are in good agreement and their correspondence 
with experimental results is investigated and discussed. 
We also study the polymer attempt frequency for crossing 
the barrier. Our studies show that the squeezing effect 
of the electric field and the lateral confinement of the 
polymer prior to its entry to the channel are very impor- 
tant factors in the height of the free-energy barrier and 
the dynamics of the polymer translocation. These effects 
were not taken into account in the previous works. We 
also study the polymer entry into an asymmetric channel 
of a-hemolysin dimensions and discuss important factors 
in the entry time and compare our results with related 
experiments. Here, in addition to extension of Ref. [l9l |. 
new results on the polymer entry into an asymmetric 
channel of the same dimensions as the a-hemolysin chan- 
nel and the polymer attempt time for passing the channel 
are presented. 

The rest of the paper is organized as follows. The 
model and the method of the simulation of polymer entry 
into an asymmetric channel is discussed in Sec. HIl Simu- 
lation results are also presented in this section. A simple 
scaling theory for statics of a polymer before its entry to 
a channel and the height of free-energy barrier against its 
entry is developed in Sec. IIIII and then compared with 
simulation results in Sec. IIVI Simulation results for the 
polymer attempt time prior to the channel entrance is 
presented and discussed in Sec. HV] Section |V] is devoted 
to the comparison of our results with related experimen- 
tal data sets published in the literature. Finally, a brief 
review of the paper and some final notes are given in Sec. 
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II. THE MODEL AND THE SIMULATION 
METHOD 

The a-hemolysin protein channel has a completely 
known crystallographic structure. It is very stable under 
the effect of the applied voltages and different chemicals. 
Thanks to such characteristics, this protein channel has 
been used in many polymer translocation experiments. 
In these experiments, it has been observed that the poly- 
mer entry time from the trans side of the channel is con- 
siderably longer than that of the cis side. For this differ- 
ence, three main reasons have been suggested: difference 
between channel vestibules size and so different entropic 
barriers against the polymer entry from the two sides, 
asymmetric distribution of charges on the channel which 
leads to polymer attraction or repulsion from the bulk 
[isj and different projections of the channel relative to 
the membrane and so entropically easier approach of the 
polymer to the cis side of the channel 17]. 

To investigate the effect of the difference between the 
size of channel vestibules, we perform MD simulations 
of electric-field-driven polymer entry into an asymmetric 
channel of a simple geometry as shown in Fig. [TJ In our 
simulations, surfaces of the flat wall containing the chan- 
nel and the channel itself are constructed by spherical 
fixed particles of diameter a. The walls are constructed 
by particles arranged on concentric circles of different 
radii and the two cylinders forming the channel (which 
are of diameters D and d) are made by particles arranged 
on consecutive circles. 3D shape of the channel in our 
simulations can be imagined by rotating the schematic 
shown in Fig. [T] around its symmetry axis. 

The polymer is modeled as a bead-spring chain. The 
beads are of diameter a and excluded volume interactions 
of the monomers with each other and with the walls and 
the channel are modeled by the shifted and truncated 
Lennard-Jones potential, 

r is the distance between the particles, e and a represent 
the strength and the length scale of the interaction, re- 
spectively and the range of interaction is rc = 2b a. The 
monomers are connected by the finite extensible nonlin- 
ear elastic (FENE) potential, 

Ubond(r) - I Q if r > i?o. 

We set the spring constant and the maximum bond 
length as fc = 70.0^ and Rq = 1.5a. 

The equations of motion are integrated with the veloc- 
ity Verlet algorithm with the step size O.OItmd, where 
TMD — f^/r" time scale. The system is 

kept at the constant temperature T = 1.0 using the 
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FIG. 2: (Color online) Polymer entry time into a cylindrical 
channel in a wall versus the channel diameter. In the first two 
points of the plot, polymer enters from one end. In the other 
points, entry in folded state is also possible. 



Langevin thermostat with the damping constant l.Or^^ 
31[. So, hydrodynamic interactions are not included in 
our model. It has been shown that consideration of the 
hydrodynamic interactions in simulation do not change 
the polymer equilibrium properties considerably. As the 
entry time is very long in the barrier-limited regime, this 
system is in a quasi-equilibrium state. 

The polymer translocation experiments are often done 
in lAI concentration of monovalent salt Q- The Debye 
screening length corresponding to such a salt concentra- 
tion is ~ O.Snm, which is smaller than the diameter of 
ssDNA, Inm. Accordingly, we assume here that there is 
no electrostatic interaction between the monomers and 
they only interact with the applied electric field. Also no 
counterions or colons are considered in our model. 

The electric field effect on the monomers is modeled 
by applying the force F = Fz to all of the monomers. 
This force is uniform in space and it is the same for all 
the monomers. All the simulations arc performed with 
ESPResSo [13. 

To measure the polymer entry time into a channel in a 
wall, we fix the middle monomer of a polymer consisting 
of 25 monomer at a close distance from the wall and leave 
the polymer to equilibrate without any applied electric 
field. Then we release the fixed monomer and turn on 
the electric field at the same time, t — Q. The polymer 
entry time, r, is defined as the time at which one of the 
monomers completely passes the channel. In these sim- 
ulations, a cylindrical constraint of diameter Dc = 14.00- 
(coaxial with the two cylinders forming the channel) is 
set around the channel to keep the polymer from diffusing 
away. The diameter of this constraint is larger than the 
radius of gyration of the polymer to avoid any confine- 
ment effects. Excluded volume interaction of monomers 
with this constraint is also given by Eq. [T] 

Presence of this constraint is justified by the fact that 
we are interested in the barrier-limited regime. In this 
regime, the time it takes the polymer to diffuse from 
the bulk to the channel is much shorter than the time 



it takes the polymer to enter into the channel. In fact, 
we assume that keeping a fluctuating polymer close to 
the pore by a constraint is equivalent to continuous dif- 
fusion of polymers of different configurations from the 
bulk to the channel vicinity. As a result, to compare en- 
try times obtained from our simulations to experimental 
ones, they should be multiplied by the ratio of diffusion 
time of the polymers from the bulk to the characteris- 
tic time of the polymer reconfiguration. One can obtain 
an estimation for this conversion factor. The character- 
istic time of polymer reconfiguration in our simulations 
is given with the Rouse time tr = ■f;f5-To [11], in which 
To is equal to the MD time scale. For a polymer with 
N — 25, this time is tr ~ tmd- For a ssDNA with 
several tens of bases, the diffusion coefficient of the poly- 
mer is D = 10~^^^. Considering the typical polymer 
concentrations in the experiments (^ fJ-M), the diffusion 
time of the polymers to the channel is around 1ms. So 
the factor for converting the simulation times to experi- 
mental ones becomes F — iiii^. 

To obtain the entry time for each given set of parame- 
ters, averages over 100-150 realizations are calculated. 
Each realization for narrower channels typically takes 
several hours with our machines. 

In Fig. [2l MD simulation result for the entry time 
of a polymer into a simple cylindrical channel in a wall 
versus the channel diameter is shown. Region 1 in this 
figure corresponds to the situation in which the polymer 
can enter the channel only by one of its ends. Region 
2 however, corresponds to the other situation in which 
entry of the polymer in folded state is also possible. Note 
that there is a two orders of magnitude difference between 
polymer entry times in regions 1 and 2 of this figure. 
Entry in the folded state is impossible when the diameter 
of the channel is so small that two monomers can not pass 
the channel simultaneously or the energy barrier for entry 
in this state is extremely high. In rest of the paper we 
only consider channels that their narrow part corresponds 
to region 1 of Fig. O 

The polymer entry time, r, from cis and trans sides of 
an asymmetric channel are shown in Fig. EJa). Chan- 
nel parameters are taken as d = a = lAa, D — 8a and 
H = 6a. In the trans side, the polymer feels no lateral 
confinement before its entry to the channel. In the cis 
side however, the polymer is confined in a wider cylinder 
before its entry to the channel narrow part. Using this 
geometry of the channel we can probe how the presence 
of a wider region before the narrow part of the channel 
affects the polymer entry time. As it can be seen in Fig. 
[31[a), existence of the wide part in the cis side reduces 
the polymer entry time relative to the trans side dramat- 
ically. 

Simulation results for a channel of dimensions very 
close to those of the a-hemolysin channel, d = 2a, 
D — 3cr, H = 5a and a = 1.5ct, are shown in Fig. [3l^b). 
Order of magnitude of the entry times and the ratio of 
entry times from the cis and the trans sides are very close 
to the experiment [l^ . It should be noted that a polymer 
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FIG. 3: (a) (Color online) Entry time, r, the time it takes 
the polymer ends to find the channel, ro, and tq exp{-^^) are 
shown respectively with the sohd, dotted and dashed lines. 
Inset: Iree-energy barrier versus the inverse of the electric 
field, (b) Entry times into a channel of the same dimensions 
as a-hemplysin. Order of magnitude of the entry times and 
ratio entry times from cis and trans sides are in agreement 
with experiment [l^. Inset: Semi- log plot of the histogram 
of the entry time, t, corresponding to data points of panel (a) 
for the trans side. Solid lines are guides to eye showing that 
semi-log plot of the histograms are almost linear. 



containing 25 monomers in our simulations is equivalent 
to a 75-nucleotide ssDNA in the experiment. For more 
details see Sec. IV Bl 

Our simulations show that the time-limiting stage in 
the process of polymer entry from the cis side is the entry 
from the wide part of the channel (the part of diameter 
D in Fig. [1]) into the narrow part. As it can be seen in 
Fig. [5J the time needed for entry in the wide part can be 
ignored. Entry time from both cis and trans sides can be 
considered as the time it takes one of the polymer ends 
to come in front of the channel narrow part by thermal 
fluctuations, tq, multiplied by inverse of the probability 
to overcome the free-energy barrier of height AC/. This 
is the well known Vant Hoff-Arrhenius law (ISl, 



Toexp 



AU 



(3) 



In situations that the polymer can enter into the chan- 
nel in folded state, tq is the characteristic time of polymer 



fluctuations in front of the channel [3J] . In our case that 
only entry of the polymer by one end is possible, tq is 
the characteristic time it takes one of the polymer ends 
to reach the channel entrance. 

In our simulations, tq and AU are measured separately. 
To calculate tq, we block the channel narrow part and 
let the polymer to equilibrate under the electric field in 
front of the channel. Then we measure the average of the 
time interval between consecutive reachings of the poly- 
mer ends to the channel. When one of the end monomers 
of the polymer falls in a range of distance from the en- 
trance of the channel narrow part, one time for reach- 
ing of the polymer ends to the channel is counted. For 
the same end of the polymer another reaching is counted 
only if it leaves the mentioned range and returns back 
again. The amount of the equilibration time depends on 
the polymer length and the electric field strength. It is 
longer for longer polymers and weaker electric fields. We 
take the equilibration time when the radius of gyration 
of the polymer gets constant and shows only small fluc- 
tuations around the equilibrium value. For every value 
of Tq, 100-150 simulations are performed. 

To measure AU, we fix one of the end monomers of 
the polymer at the entry of the channel and let the rest 
of the monomers to equilibrate under the electric field. 
Then we release the end monomer and probe whether 
it enters the channel narrow part or returns back. The 
logarithm of the probability of successful entry (ratio of 
the number of successful entries to the total number of 
simulations) equals to Number of realizations to 

obtain this quantity is typically lO'^-lO^ in our simula- 
tions, depending on the barrier height. 

As it is shown in Fig. El^a), values of both tq and 
AU for the cis side are smaller than those of the trans 
side. In Fig. [Sja), product of tq and exp ^^^^ is also 

compared with r, which shows that they are of the same 
order of magnitude. In this figure, the size of error bars 
on r and tq data points are of the same order of the data 
itself. Such large error bars are resulted from the expo- 
nential distribution of these quantities (see the inset of 
Fig. EJb)). For any variable of exponential distribution, 
it is known that the mean and the standard deviation 
are of the same value. Actually, exponential distribution 
of entry time originates from the fact that consecutive 
pol ymer entries into the channel are independent events 



III. THEORY 

In this section, we present a simple scaling theory 
for the model system introduced in the previous section 
which will be tested with the simulation results in the 
next section. An important point in this theory is that 
the electric field compresses the polymer to the wall be- 
fore its entry to the channel. Investigation of the polymer 
conformation in this compressed state is necessary for 
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FIG. 4: Schematic of a polymer under the electric field of 
strength Eo behind a channel of diameter d in a wall, (a) 
The polymer can be viewed as a 2D chain of blobs of size Lj . 
(b) A polymer which is also confined laterally by a cylinder 
of diameter D. The polymer can be considered as a chain of 
blobs closely packed inside a cylinder of height L//. (c) A 
segment of the polymer entered to the channel inside which 
the electric field is of strength Ei. 



obtaining the free-energy barrier, AU, and the polymer 
attempt frequency, Tq^. 

As case /, consider a polymer with N monomers of 
length b and charge q, which is pushed to a wall by an 
applied electric field. The electric field tends to move 
the polymer toward the wall and decrease its electric en- 
ergy. On the other hand, entropy resists the polymer con- 
finement in a small volume near the wall. Competition 
between these two effects determines the polymer layer 
thickness, Lj (see Fig. |4](a)). Actually, the electric field 
effect on the polymer is similar to the polymer adsorption 
to the wall, which is modeled with the polymer confine- 
ment between two walls [35i] . Here, one can introduce the 
blob size, a length scale in the system smaller than which 
the polymer has its unperturbed statistics and beyond it, 
the confinement by the wall dominates. So, the polymer 
can be modeled as a 2D chain of blobs of size Lj. Inside 
these blobs the Flory statistics reads Lj ~ bgj in which, 
gi is the number of monomers in each blob. Free-energy 
of the polymer under the electric field is of the order of 

keT per blob; F.^t ~ fesTf - ksTN (;^) " . 

As case //, consider a similar polymer which is also 
confined from sides by a cylinder of diameter D (see Fig. 
|4](b)) . In this case also, the effect of the electric field is 
equivalent to confining the polymer between two parallel 
walls. Because of the simultaneous confining effects of 
the electric field and the cylinder, this situation is similar 
to a polymer confined in a closed cavity. In this case, 
blobs can be defined that inside them the Flory statistic 
governs, ^ ^ bgjj. These blobs are closely packed inside 
the cavity, y^'^ ~ fi. Here ^ and gn are the blob size 
and the number of monomers in each blob, respectively, 
and fl ^ LjjD^ is the cavity volume [sB]. These two 
relations give the number of monomers in each blob gn ^ 

(/fp')^""^' ^^'^ ^^"^ polymer free-energy is obtained as 



We obtain the electric energy of the polymer layer, 
feicj in both cases with a mean field approximation. 
A uniform distribution is assumed for monomers in the 
layer of thickness {m = I, II)- Taking zero of the 
electric energy on the wall, we obtain Feic ~ NqEoLm- 

Because of the long entry time, polymer reaches to an 
equilibrium state before its entry to the channel. So we 
can find the thickness of the polymer layer in both cases 
by minimizing the sum of electric and entropic energies 
of the polymer, Fent + Feic, with respect to Lm- 

„..(^)*,._.(g^)^(^)*, 

(4) 

Free-energy per monomer before entry to the channel 
is ksT / gm and the number of monomers in each blob 
is obtained as 

Consider a polymer in each one of cases / and // that 
enters a channel of diameter d on the wall. Inside the 
channel, we take the polymer as a chain of blobs of diam- 
eter d (see Fig. 111(c)). Statistics of the polymer segment 
inside a blob is unperturbed and d ^ bg^^ in which g^ is 
the number of monomers in each blob inside the channel. 
Entropic energy of each monomer inside the channel is 

^ ■^^j where g^ ^ (f ) • A chain of n monomers has 
a length I ^ —d n(^)^/'^(i inside the channel. 

Hence, total change in the polymer free-energy af- 
ter entry of n monomers into the channel is Ai^„i '~ 

kgT {{i^y - j-^n - n^qEii^y/^d. In this relation, 

the term linear in n comes from the entropic energy and 
the quadratic term is the electric energy. Maximum of 
AFjn with respect to n gives the height of the free-energy 
barrier, AUm, which determines the rate of the polymer 
entry into the channel: 

AUm ( ^(dV^" 

keT ^ [\dj gm) qE,d \b ) ' ^ ' 

The main dependence of AUm on E comes from the 
E~^ term. The prefactor which contains g has a weak de- 
pendence on the electric field and over small intervals of 
the field strength this term is almost constant. Such lin- 
ear behavior of the energy barrier with E~^ can be seen 
in the inset of Fig. EJa). But over large electric field in- 
tervals, the changes coming from g-dependent prefactor is 
considerable and AU versus E~^ would show deviations 
from the linear behavior. This has also been observed in 
the experiment [st} . 

In the inset of Fig. [3Ka) , the slope of the fitted line for 
case / is larger than that of case //. This is because of 
the g-dependent prefactor in Eq. [5] and larger g for case 
/. 
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Equation |6] can also be applied to free-energy barrier 
in entropic trapping experiments (37j . In these experi- 
ments, the mobility of the charged polymer under electric 
field is measured in a channel with alternative wide and 
narrow regions. This mobility depends on the polymer 
trap time at the entrance of the narrow regions. Trap 
times are determined by the free-energy barrier against 
the entry of the polymer from the wide parts into the 
narrow parts of the channel. This method has been used 
for polymer separation 38]. It has been shown that in 
the experiment conditions, the free-energy barrier does 
not depend on the polymer length, iV, and dependence 
of the trap time on N comes only from the dependence of 
the polymer attempt time for passing the barrier on the 
polymer length. This time decreases with increasing the 
polymer length [2l|, [s^, . The interesting point in Eq. 
|6] is dependence of the free-energy barrier on the poly- 
mer length in case //. In this case, AU decreases with 
increasing the polymer length, N. Thus by tuning the 
channel parameters the situation for the polymer before 
its entry to the channel narrow part could be similar to 
case //. In this way, dependence of the trap time on the 
polymer length would get stronger and may lead to more 
optimized polymer separation. 

The number of monomers in each blob, gm, can not 
be smaller than 1 or larger than the total number of 
the monomers, N. In weaker fields (or shorter polymer 
lengths), the field is not strong enough to confine the 
polymer in a layer close to the wall and gm > N. As a 
result, the polymer takes its unperturbed conformation 
and the radius of gyration is i? ^ bN"^ . In the other limit, 
g < 1 means that the polymer has lost all of its entropy 
under the electric field and it enters the channel without 
feeling a free-energy barrier. These two limits have also 
been observed in experiments. A crossover voltage below 
which no entry happens and another voltage limit above 
which the entry time chan ges linearly with the applied 
voltage have been reported (23l [27j . 



IV. SIMULATION RESULTS 



A. Statics of the polymer before entry to the 
channel 



In simulations for checking the theoretical results on 
the polymer statics before its entry to the channel in 
both cases / and //, we block the channel entrance and 
let the polymer to equilibrate under the electric field. 
Then we measure the radius of gyration of the polymer 
parallel and perpendicular to the wall and find their scal- 
ing exponents with the polymer length and the electric 
field strength. 

In case /, the polymer feels no lateral confinement on 
the wall. In this case, the blobs can be considered as the 
effective monomers of a 2D polymer. The blobs have a 
self-avoiding random walk on the wall and the radius of 
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FIG. 5: (Color online) (a) Exponents /3 and 7 in relations 
L oc A'^'^ and R\\L^ oc A'^^ for three different values of the field 
strengths, (b) Exponent q in L oc E'" for three different 
polymer lengths. 



gyration of the polymer parallel to the wall is 



HI 



(7) 



In this relation, gj is substituted from Eq. [5] and | is the 
Flory exponent in 2D. If the diameter of the confining 
cylinder is larger than this radius of gyration, Z) > , it 
has no confining effect on the polymer and the situation 
is similar to case /. But if the cross section of the cylin- 
der is smaller than the monomers total area, < Nb'^, 
the monomers pile on each other and case // is relevant. 
Thus, for a given diameter of the confining cylinder, short 
and long polymers (or weak and strong fields) are equiv- 
alent to cases / and //, respectively. One expects the 
crossover between these regimes to occur in the interval 



^)^<iV<(f)^when 



qEb 



1. 



o / \ o / ' R _t 

According to Ec(s. |4]and[7l one can write the relations 
_L oc and R\\L^ oc N'' for dependence of perpendicu- 
lar and parallel components of the radius of gyration of 
the polymer on the number of its monomers. Then, it 
is obtained that /?/ = and 7/ = 0.75 for case /, and 
= ~ 0.6 and 7// = ~ 0.14 for case //. These 
exponents which are calculated from the local slopes of 
the log-log plots of L and R\\Li versus N are shown in 

Fig. [HJa). For the electric fields = h^TE ^^^'^ 
creasing N, crossover from case / to case // can be seen. 
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The weaker electric field = ^ is out of the validity 
liniits of the theory and we have 7~/3~j/ = 0.6 for all 
the polymer lengths. In these simulations, the diameter 
of the confining cylinder is 13 = 12a. With this diame- 
ter, shorter polymers do not feel the confining effect of 
the cylinder (case /). Longer polymers however, corre- 
spond to case //. The crossover between the two cases 
happens at 27 < < 144, in agreement with the theory. 

Using Eq. m dependence of the radius of gyration of 
the polymer perpendicular to the wall on the electric field 
strength can be written as i oc i?^". For cases / and //, 
ai — 0.37 and an — 0.44 are obtained. This exponent is 
calculated from the local slope of the log- log plot of L ver- 
sus E (see Fig. [SJb)). For the polymer length N — 256 in 
the weaker fields, the exponent a is very close to the ex- 
pected value 0.4 (note that values of the exponent a are 
very close to each other in these cases and distinguishing 
between them is not possible in the simulation). With 
increasing the electric field, the exponent a decreases, 
which is the result of approaching the validity limit of 
the theory, g < I. Assuming that the scaling constant in 

Eq. Elis of the order 1, gn = = 1 gives the 

0.5 for this limit. For the two shorter poly- 



value 

Kb 1 



mers, it is seen that in the range of 0.03 < < 1 fo'" 
iV = 64 and in the range of 0.3 < §^ < 1 for = 16, 
the exponent a is between 0.3 and 0.5 and the curve slope 
is smaller. But despite the theory, no plateau can be seen 
in this curve, which may be a result of the failure of the 
mean field approximation. 



B. Free-energy barrier 




NqEb/kgT 



140 



FIG. 6: (Color online) (a) , (b) Symbols are simulation results 
for free-energy barrier versus the electric field for different 
values of the polymer length in cases / and II. In (b) the 
axes are scaled such that all the data sets fall on a master 
curve. The solid line is the fitted function of the theory, Eq. 
[6l to the simulation data. Insets show free-energy barrier 
versus . Deviations from linear behavior can be seen over 
large intervals of the electric field, specially in case II. 



To check the main result of our scaling theory, Eq. |6l 
we calculate the height of the free-energy barrier with the 
same method as described in Sec. HIl for different values 
of the polymer length and the electric field strength (see 
Fig. ini). The diameter and the length of the channel 
in these simulations are d = 1.4(7 and H — 6a and the 
diameter of the confining cylinder is D = 8a. As it can 
be seen in this figure, our scaling theory can explain the 
exponents but not the coefficients. We use the function 
(a - bx^-^)^/x for fitting to ^ versus |^ for case / 

and the function {a — cx'^'^)'^/x for fitting to ^f^j, versus 

for case //. a, b and c are fit parameters. The 
parameter a depends on the ratio of the monomer size 
to the diameter of the channel. It is taken the same 
in separate fittings to the data sets of cases / and //. 
As it can be seen, there is a good agreement between 
simulation and theory and the fit parameters are close 
to their expected values. Taking d &, the expected 

values are a - (^)^^^ - 1, 6 - (^)"^^^ 1 and 

(s) ^" {2) ^ ^ ^'^ reasonably well agreement 
with the values obtained from fitting to the simulation 
data, a = 3.02, b = 1.50 and c = 0.32. 



In the insets of Fig. [5] (a) and (b), deviation of AU 
versus E^^ from linear form in large intervals of -E"^, 
no dependence of the free-energy barrier on the polymer 
length in case I and its decrease with the polymer length 
in case // can be seen. 



C. Polymer attempt time for crossing the barrier, 

To 

According to the Rouse model, internal monomers of 
the polymer have an anomalous diffusion and their mean 
square displacement changes with time as (Ar^) oc Af°-^ 
[35|. In our model, one can consider the end monomer 
of the polymer as a particle that is anomalously diffus- 
ing in a closed cavity of volume 57 = LD^, with L the 
radius of gyration of the polymer perpendicular to the 
wall. The mean polymer attempt time for passing the 
barrier is then equal to the mean first passage time for 
this particle. This time changes with the 4th power of 
the system characteristic length, Qs; tq oc fis [40|. Our 
simulation results for the attempt time are very limited 
and scattered and we cannot reject or prove this relation. 
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FIG. 7; (Color online) Distribution of z component for one 
of end monomers of the polymer for cases I and II. The 
values of z components corresponding to the maximum of the 
distributions are Zmj = 1.55 and Zm,ii = 1.15. Inset: tq 
versus f2 = LD^ . Although the values of for cases / and 
II are close to each other, the values of tq for the two cases 
are very different. It may be resulted from the failure of the 
mean-field assumption for the distribution of the monomers 
in our theory. 

But a point which can be emphasized here is the large 
difference between the polymer attempt times from the 
cis and trans sides, despite the similar values of 51 for the 
two sides (see the inset of Fig. [T]). Actually, in the trans 
side, the end monomers do not get close to the wall eas- 
ily for entropic reasons. But in the cis side, the polymer 
is confined in a closed volume and the pressure of the 
other monomers push the end monomers toward the wall 
as is shown in Fig. [71 In this figure, there is a pick in 
distribution of the position of the end monomer on the 
wall in case //, which is not seen for case /. This point 
can explain the large difference between the two cases 
and shows the failure of the mean field approximation 
for this part of the problem. 

V. COMPARISON WITH EXPERIMENT 
A. Partition coefficient 

In the experiments, the partition coefficient is mea- 
sured to calculate the confinement free-energy of the 
polymer. The partition coefficient is the ratio of polymer 
concentrations inside the channel and in the bulk solu- 
tion. A voltage difference is applied to keep the polymer 
in the channel and ion current through the channel is 
measured [U 14^. Also, quantities such as polymer re- 
action rate with a residue inside the protein channel can 
be measured Q. 

From our theory, we expect the confinement free- 
energy of a polymer of length N inside a narrow channel 
of diameter D under the electric field E to obey the re- 
lation ^ - ^ cx N^-^E°-^ (■I.V'^. This relation is 
valid for high salt concentrations (no long range electro- 



FIG. 8: (Color online) Logarithm of capture time versus AV 
extracted from referenced articles (symbols) and fitted equa- 
tion of our theory, Eq. [S] (solid lines). All the data sets 
are from polymer translocation experiments through the a- 
hemolysin channel. As it can be seen, all data sets are followed 
well by the suggested function. Fit parameters are shown in 
Table [TTl 

static interaction between the monomers) and for poly- 
mer lengths and electric field strengths that are in the 
validity range of the theory. 

There are numerous suggestions for dependence of the 
confinement free-e nerg y on the polymer length in the lit- 
erature. In Ref. [4l| exponents 1.6, 1.35 and 3.2 are 
obtained for confinement inside different channels. Ex- 
ponent 3.2 is explained in Ref. [36j and it is attributed to 
the almost closed shape of the channel used in the exper- 
iment. The two other exponents are in good agreement 
with our theory. In Ref. Q it has been shown that the 
partition coefficient obeys the de Gennes theory. In this 
theory, the confinement free-energy of the polymer inside 

a cylinder without any electric field is ^ ^ {.Ts)^ ■ 
The linear dependence of the free-energy on the poly- 
mer length in this experiment may be because of the low 
salt concentration, lower applied electric field or smaller 
polymer length that are beyond the polymer compressing 
limit by the electric field and hence beyond the validity 
range of our theory. 

B. Dependence of the entry time on the applied 
voltage 

According to Eq. [31 the entry time depends on the free- 
energy barrier height as well as the polymer attempt time 
for passing the barrier. Considering the weak dependence 
of To on the applied voltage (see Fig. [3lja)), dependence 
of the entry time on the applied electric field for both 
cases / and II can be written using Eqs. [3l and [Bias 

T], A and C, are parameters independent of the applied 
voltage. Most of the applied voltage drops along the 
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TABLE I: Parameters of the three polymers used in the ex- 
periments and the corresponding course-grained parameters. 
The polymers are single-stranded DNA, dextran sulfate and 
poly(styrene) sulfonate. 





ssDNA 


DS 


PSS 


real monomers separation (A) 


3.3 


2.5 


2.5 


real monomers charge (e) 


1 


2 


1 


polymer diameter (A) 


10 


7 


8 


coarse-grained monomer size (b) (A) 


10 


7 


8 


coarse-grained monomer charge (q) (e) 


3 


5.6 


3.2 


polymer persistence length (A) 


15 


16 


14 



channel and dependence of the electric field inside the 
channel on the voltage reads Ei — where H is the 
channel length. Here, it is assumed that the electric field 
outside the channel is uniform and has a linear relation 
with the field strength inside the channel, Eo — SEi, 
where i5 is a constant parameter. We obtain b and q in Eq. 
[8] from coarse-graining rule of the polymer. We take the 
coarse-grained monomers size, 6, the same as the polymer 
diameter in the experiment. The charge of coarse-grained 
monomers, q, is determined from known linear charge 
density of the polymer. Parameters related to several 
flexible polymers and their coarse-grained monomers are 
shown in Table HI 

According to Eqs. [3] and H C Inro, - (f)^'^^ 
for both cases / and //, A ^ (|) ^ ^" 5^-^ for case / 

and A ~ (f )^~^ <5°-^ (^)°^ for case //. We use these 
quantities as fit parameters to the experimental results as 
shown in Fig. |51 All these results are related to polymer 
passage through the a-hemolysin channel, so H = lOnm. 
Fit parameters are shown in Table |lll 

Because of the complex shape of the a-hemolysin chan- 
nel, the profile of the electric field strength in the chan- 
nel vicinity is not very clear in these experiments and we 
can only investigate the fit parameters qualitatively. C 
depends on tq and hence on the polymer concentration 
in the solution and the polymer length, r/ depends on 
the ratio of diameters of the polymer and the channel. 
In addition, 77 could be dependent on the electrostatic 
interactions between the polymer and the channel. Re- 
garding almost the same charge and diameter of the poly- 
mers used in the experiments listed in Table HIl and also 
considering the same protein channel used in these exper- 
iments, the same value of rj for all data sets is acceptable. 
A is the main parameter of interest here. It is a measure 
of polymer compression before its entry to the channel. 
Diameter of the a-hemolysin channel in the trans side is 
smaller than the cis side. This causes a higher resistance 
for this part of the channel and a higher voltage drop 
and stronger electric field inside it. The stronger electric 
field leads to stronger compression of the polymer. It 
can be seen that A is larger for the data sets which are 
for translocation from the trans side compared to those 



TABLE II: Parameters (, rj, and A from fitting of Eq. [8] to 
experimental data. Result is shown in Fig. [5] 





Experimental data 




V 


\ 

A 


1 


ssDNA, cis side (Ref. [23]) 


-5.42 


3.14 


0.60 


2 


ssDNA, cis side (Ref. [18]) 


-4.33 


2.38 





3 


ssDNA, trans side (Ref. fl8]) 


0.27 


3.22 


1.77 


4 


ssDNA, cis side (Ref. [26]) 


-3.01 


3.11 


0.90 


5 


DS, cis side (Ref. [27]) 


-7.09 


2.35 


0.01 


6 


DS, trans side (Ref. \2^) 


-6.13 


3.85 


0.51 


7 


PSS, cis side (Ref. [29|) 


-6.01 


2.83 






for the cis side. In addition, in translocation from the 
cis side, when the applied voltage is higher, compression 
effect is stronger and the value of A is larger. The data 
sets 1 and 4 are for high voltages and data sets 2, 5 and 7 
are for low voltages (all of them are for the cis side) . The 
point which should be noted here is that according to 
Ref. [2^, in translocation from the cis side, the polymer 
is not compressed inside the wide part of the channel. 
So, entries from both cis and trans sides may correspond 
to case / of our theory, not case //. 

VI. SUMMARY AND DISCUSSION 

Entry of a flexible polymer into a narrow channel which 
is the most time-consuming stage of the electric-field- 
driven polymer translocation through nanoscale channels 
in the barrier-limited regime has been studied theoreti- 
cally and using MD simulations. Two different geome- 
tries for the channel entrance have been considered to 
investigate the effect of the lateral confinement of the 
polymer before its entry to the channel narrow part. The 
height of the free-energy barrier has been obtained as a 
function of the strength of applied electric field and the 
length of the polymer for the two entrance geometries. 
There is a very good agreement between our theoretical 
and simulation results. Our suggestion for dependence of 
the polymer entry time on the applied voltage across the 
channel fits with the related experimental data and the 
behavior of the fit parameters are described. Our theory 
also explains the polymer partitioning experiments. 

With the aim of describing the reported difference be- 
tween polymer entry times from cis and trans sides of 
the asymmetric a-hemolysin channel, MD simulation of 
polymer entry into a channel of dimensions very close 
to those of a-hemolysin has been performed. The order 
of magnitude and the ratio of polymer entry times from 
cis and trans sides obtained from our simulations are in 
agreement with the corresponding experiment. Accord- 
ing to our results, existence of a confining space prior 
to the channel narrow part in the cis side, very similar 
to the geometry of the a-hemolysin channel, causes the 
same difference between the entry times from the cis and 
the trans sides as has been observed in the experiments. 
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Free-energy landscape seen by the polymer when it enters 
from the cis side is shown schematically in Fig. [T] When 
the electric field is strong enough, after passing one step 
of the energy barrier (entering into the wide part), the 
polymer is trapped in a region where its free-energy is 
lower relative to the outside of the channel (the gray re- 
gion in the lower panel of Fig. [IJ . It causes the passage 
from the two steps of the energy barrier to happen as 
two separate stages and the passage probability from the 
whole barrier not to be the product of passage probabili- 
ties from the two individual steps. If the diameter and/or 
the length of the wide part or the electric field strength 
is small enough, the region in which the polymer has a 
lower free-energy disappears. Then the total height of 
the free-energy barrier would determine the polymer en- 
try time. We have observed such behaviors in our test 
simulations. 

The channel projection out of the membrane has been 
shown that has not an important role in the polymer 
entry process from the cis side in our simulations. In 
our test simulations with the real structure of the a- 
hemolysin channel taken from the protein data bank (43| , 
we observed that the projected surface of the channel in 
the cis side has a larger diameter than the radius of gy- 
ration of the polymer (for a polymer of length N — 25). 
So the projected part of the channel effectively interacts 
as a wall with the polymer. Therefore, in all our simu- 
lations we have put the wall at the channel ends in both 
the cis and trans sides. In addition, in Ref. flT'l that the 
channel projection is suggested as a reason for difference 
between entry times from the cis and trans sides, the 
electric field outside the channel is taken to be zero. In 
our simulations however, electric field outside the chan- 
nel compresses the polymer and the entropic cost of the 



attachment of the polymer end to the wall is small. The 
effect of the projected part may be important at weaker 
electric fields or for longer polymers. 

To make the model used in our theory and simula- 
tion more similar to the real system in the experiment, 
two points should be considered. First, the electric field 
outside the channel should be calculated from more ac- 
curate theories jl^ and also the field gradient inside the 
a- hemolysin channel should be taken into account p3 |. 
Second, some percentage of the polymers that enter from 
the cis side of the a-hemolysin channel are compressed 
inside the wide part of the channel before their entry to 
the narrow part j23l |. In this situation, a portion of a 
long polymer would remain out of the wide part of the 
channel. This case is not studied in our theory. In this 
situation, the polymer is divided into two parts: one is 
compressed inside the cylinder and the other is free on 
the wall. The part outside the channel may act like a 
spring which drags the inner part out of the channel |45| . 

In many of the polymer translocation experiments, 
double stranded DNA which its bending rigidity is con- 
siderably larger than ssDNA has been used. A more com- 
plete theory is needed to be developed for semiflexible 
polymers. More complex regimes have been predicted 
for semifiexible polymers confined in a nanoscale space 
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